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We report on an astonishing switching synchronization phenomenon in one-dimensional mem¬ 
ristive networks, which occurs when several memristive systems with different switching constants 
are switched from the high to low resistance state. Our numerical simulations show that such a 
collective behavior is especially pronounced when the applied voltage slightly exceeds the combined 
threshold voltage of memristive systems. Moreover, a finite increase in the network switching time 
is found compared to the average switching time of individual systems. An analytical model is 
presented to explain our observations. Using this model, we have derived asymptotic expressions for 
memory resistances at short and long times, which are in excellent agreement with our numerical 
calculations. 


PACS numbers: 64.60.aq, 73.50.Fq, 73.63.-b, 84.35.+i 

Introduction - Synchronization is the term that is fre¬ 
quently used to describe the coherent dynamics of an 
ensemble of interconnected dynamical units, namely, dy¬ 
namical units forming networks. The networks are ubiq¬ 
uitous in nature and technology, and, therefore, it is not 
surprising that the phenomenon of synchronization has 
been studied and observed in a wide range of dynamical 
systems. These systems include, for example, neurons 
BE], power grids [ 3111 ], coupled lasers [5|, oscillators, 
social systems |5], etc. A lot of attention has be drawn 
to the synchronization of chaotic systems mis] - an in¬ 
triguing emergence of collective dynamics of a number of 
chaotic units linked with a common signal or signals. Os¬ 
cillator networks |9] are another example of widely stud¬ 
ied systems with synchronization. 

In this Letter, our attention is focused on memris¬ 
tive (memory resistive) networks. These networks are 
composed of individual memristive elements uni, which 
now are of considerable interest for a variety of applica¬ 
tions. In these passive resistive electronic devices, the 
resistance depends on the history of signals applied. The 
ability of memristive systems to store and process infor¬ 
mation on the same physical platform makes them ideal 
for unconventional computing applications dUlIl]. In 
fact, boolean logic operations with small memristive net¬ 
works were experimentally demonstrated few years ago 
m- Moreover, it was theoretically shown that larger 
memristive networks could solve maze m and shortest 
path optimization m problems in a single step compared 
to multi-step algorithms employed in conventional com¬ 
puters. Therefore, it’s of a real importance to understand 
the dynamical properties of memristive networks. 

Recently, two of us (VAS and YVP) have found that in 
one-dimensional memristive networks subjected to adia- 
batically increasing voltage, the effective switching rates 
of memristive systems strongly depend on their polari¬ 
ties [T6|. It has been demonstrated (on the level of in¬ 
dividual memristive elements) that an abrupt (acceler¬ 


ated) switching occurs when the memristance (memory 
resistance) of a given memristive system in the network 
increases at the given voltage polarity. A slow (decel¬ 
erated) switching takes place in the opposite case [l6] . 
However, this prior work leaves open the question of the 
switching behavior beyond the adiabatic limit, namely, 
when the applied voltage is initially high enough to in¬ 
duce the dynamics of several memristive systems. This 
is precisely the aim of the present Letter, which explores 
the switching dynamics of one-dimensional memristive 
networks subjected to sufficiently high voltages. Accord¬ 
ing to our findings, an interesting switching synchroniza¬ 
tion effect takes place when all memristive systems switch 
from the high to low resistance state. Our consideration 
of the switching synchronization effect employs both nu¬ 
merical and analytical techniques. 

Fig. [2 shows a one-dimensional memristive network 
connected to a constant voltage source V. It is as¬ 
sumed that the polarities of all memristive systems are 
the same, and, at the initial moment of time t = 0, all 
memristive systems are in the same high resistance state 
Ri{t = 0) = Roff- Such configuration is similar to the 
initial configurations in the maze |14j and shortest path 
problem m solving networks and thus of importance for 
potential applications of memristive systems. 
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FIG. 1: (Color online) One-dimensional network of N mem¬ 
ristive systems Mi connected to a voltage source V. Here we 
assume that all the memristive systems Mi are connected with 
the same polarity and Ri{t = 0) = Roff, where Roff is the 
high resistance state. 
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We assume that Fig. network employs voltage- 
controlled memristive systems with threshold that have 
been experimentally realized with different materials 
combinations m In our numerical simulations and an¬ 
alytical calculations presented below, we use a model 
of voltage-controlled memristive systems with threshold 

mm 

I = R-\x,)V. (1) 

dxi 

dt 

where / and Vi are the current through and the volt¬ 
age across i-th memristive system, respectively, Xi is the 
internal state variable playing the role of memristance, 
Ri(xi) = Xi, /3i is a, positive switching constant charac¬ 
terizing the intrinsic rate of memristance change when 
\Vi\ >Vt, Vt is the threshold voltage, and -I- or — sign is 
selected according to the device connection polarity. Ad¬ 
ditionally, it is assumed that the memristance is limited 
to the interval [Ron, Raff] (note that Ron < Raff)- 

Numerical Results.- Let us, first of all, consider switch¬ 
ing in an ensemble of individual memristive systems char¬ 
acterized by a probabilistic distribution of the parameter 
Pi- For the sake of simplicity, we consider a flat random 
distribution of this parameter keeping all other param¬ 
eters of memristive systems the same. Fig. ia) shows 
the time-dependence of memristances in an ensemble of 
memristive systems each subjected to the same voltage 
Vi = 1.05Vt, i = 1,.., N, which is the average voltage per 
system in the network considered in the next paragraph. 
It’s not surpising that the switchings of these memristive 
systems occur at very different rates defined by specific 
individual values of switching constants Pi. 

Next, we consider the collective switching, namely, 
the switching of memristive systems that form one¬ 
dimensional networks as the one sketched in Fig. [l] 
Fig. [Kb) shows the main result of this Letter. This 
plot demonstrates a surprising switching synchronization 
effect in which the effective switching rates of unlike sys¬ 
tems become the same. This is a truly remarkable be¬ 
havior that has not been anticipated in the literature. 
The basic principles of this behavior are related to the 
phenomenon of the decelerated switching m that, un¬ 
like our previous investigation [IB]. takes place simulta¬ 
neously in every component of the network. Technically 
speaking, the switching of memristive systems with larger 
values of Pi can not proceed fast as the decreases of their 
memristances also suppress the voltage falls across them. 
At the same time, the voltage falls across memristive sys¬ 
tems with smaller values of Pi increase compensating the 
smallness of their Pi. As a result, the switching of all 
memristive system occurs coherently with approximately 
the same effective rate. 

Fig. Kb) also demonstrates that the total switching 
time for the entire network is longer than the switching 


j±sign(I4)/3i(|Fi| - Vt) if \Vi\ > Vt 


1» 


otherwise 


( 2 ) 


T-'-1-'-r 



0 1000 2000 3000 4000 5000 

Time (in units of 



Time (in units of r^) 


FIG. 2: (Color online) The memristance Ri{t) of A = 30 
memristive systems with Roff/Ron = 100 (a) individually 
subjected to the same voltage Vi — l.OSFt, and (b) forming 
ID network (as in Fig. Q subjected to F = l.OSAVt- These 
plots have been obtained with a random flat distribution of 
parameters Pi in the interval [(/3) — ^Pl‘^, {P) + A/3/2]. The 
time is measured in units of tq = Ron/i{P)Vt). 


time defined by the average value of Pi-s, which is (/3). 
This general type of behavior has been verified for many 
different realizations of memristive systems. Fig. Kshows 
some additional data points (extracted from numerical 
simulations) demonstrating a monotonic increase in the 
network switching time with distribution width. In fact, 
according to our analytical theory presented below, the 
network switching time is actually proportional to (1//3) 
instead of 1/{P) (See Eq. [H]). For a flat random distri¬ 
bution of Pi in the range [(/3) — NPI2, (/3) -h A/3/2] and 
iV» I 


I 


(/3) + ^ 


(1//3) = ^In 

A/3 (p)- 
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2 


( 3 ) 


The dashed curve in Fig. [^ shows the perfect agreement 
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FIG. 3: (Color online) Network switching time as a function 
the distribution width A/3. The switching times have been 
calculated for a network of N = 1000 memristive systems with 
a flat random distribution of switching constants subjected to 
V — l.OSNVt. The dashed curve is plotted assuming that the 
switching time is proportional to (1//3) given by Eq. (|^. 


of Eq. <© with our numerical results. 

If the distribution of Pi-s is not flat then one can show 
that for any distribution of Pi the difference {l/P) — 
1/(|S) > 0. Additionally, this difference normally grows 
with the distribution width. For example, if all odd cen¬ 
tral momenta are negative or equal to zero (as, for exam¬ 
ple, for the Gaussian distribution), then this difference 
cannot be less than ((/3^) — (/?)^)/(/3)^. Thus, our obser¬ 
vation of the switching time increase is valid on average 
for any distribution of Pi. 

Analytical Model. — If the initial memristances are 
the same and the applied voltage exceeds the combined 
threshold voltage NVt , then one can realize that the volt¬ 
age fall across any memristive system exceeds its thresh¬ 
old voltage at any time. Moreover, in the case of a dis¬ 
tribution of initial memristances, the same is true either 
from t = 0 or after an initial equilibration period. There¬ 
fore, in the region of parameters of interest, Eq. (§ can 
be generally written as 


RPt) = -P,[Vpt)-Vt], (4) 

where i = 1,...,N, Vpt) = VRi{t)/R{t), and R{t) = 

N 

^ Ri (t) is the total memristance. The minus sign in Eq. 

i—1 

0 takes into account the selected connection polarity 
of memristive systems in the network (such that their 
memristances decrease at positive applied voltages). 

Let us search for the solution of Eq. 0 in the form 

RPt) = Cpt)e , (5) 

where C'i (t) is a time-dependent function and the integral 



FIG. 4: (Color online) Comparison of asymptotic expressions 
(Eqs. (|^ and Q) for the time-dependence of the total mem¬ 
ristance R{t) and numerical (exact) solution. The numerical 
solution (solid black curve) has been obtained for the same 
realization of memristive systems, model and simulation pa¬ 
rameters as in Fig. Eq. ([^ curve is plotted in the linear 
approximation. 
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in the exponent is actually the charge (/ E/i?“^(r)dr = 

0 

q{t)) flown through the network by the time t. Substi¬ 
tuting Eq. 0 into Eq. 0 > one can find the following 
expression for Ri{t): 


RPt) = i?,(0)e-^‘«(‘) + J (6) 

0 

Taking into account that i?i(0) = Roff, the sum of Eqs. 
^ yields 



As q{t) can be expressed through R{t) (see the definition 
of q{t) below Eq. 0), Eq. (j^ can be considered as the 
nonlinear integral equation for R{t). 

While it is difficult to find the exact solution R{t) from 
Eq. 0. this equation can be effectively used to derive 
the asymptotic behavior of R{t) in the most important 
limiting cases. In particular, in the short time limit, one 
can expand R{t) = RQ{l — at+bt‘^+0{t^ j) using unknown 
constants a and b, and get q{t) = V/2+0{t^))/R q. 
Using Eq. Q one can find 

R(t) = NRoff-{P)SVt+ ^fJ^^^" +0(f^), t ^ 0, (8) 

2A Roff 

where SV = V — NVt is the voltage excess above the com¬ 
bined threshold, and DP = (/3^) — (/3)^ is the dispersion 
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of switching constants Pi. While the first and the second 
terms in Eq. (|^ correspond to the averaged effect of evo¬ 
lution of individual memristive systems, the third term, 
being proportional to the dispersion, is always positive 
and describes effect of collective evolution of memristive 
network. Note that the expression ([^ is valid only when 
the second and third terms are small compared to the 
first one. 

A different asymptotic expression can be found in the 
long time limit, namely, when Piq(t) 3> 1. This limit also 
implies the optimal synchronization condition 5V <C V 
as demonstrated below. When Piq{t) ^ 1, the main 
contribution to the right-hand side of Eq. Q comes from 
the upper limit of the integral with respect to P. Using 
this observation one can derive the following main term 
of the long time asymptotic 


R{t) = {NRoff-PHSVt){l + 0{SV/V)), SV^+0, (9) 
where Ph = (l//3)“^. 

To specify the applicability conditions of Eq. 0 , one 
can calculate q{t) using Eq. Then the condition 

Piq{t) 3> 1 can be presented as 


P^V 

PhSV 


In 


NR. 


off 


m 


> 1 . 


( 10 ) 


Eq. (10) can be sub-divided into the optimal synchro¬ 
nization condition dV V (also observed in our numer¬ 
ical studies) and the condition of long times such that 
the total resistance R{t) is much less than its initial value 
NRoff. The total switching time T for the network can 
be easily computed substituting R{T) = NRon in Eq. 
0. This gives 

^i^off — Ron) / 1 


T = 


V-NVt 


-^hsv«v. 


(11) 


Eurthermore, in the typical situations when Roff ^ Rom 
Ron in the nominator of Eq. 0 can be omitted. 

Fig. shows a comparison of the numerically obtained 
solution with the asymptotic expressions given by Eqs. 
([^ and 0. Clearly, the asymptotic expressions are in 
the excellent agreement with the numerical solution for 
R{t). 

It is interesting to note that Eq. 0 also delivers a good 
approximation for all times when SV V. This allows to 
find the approximate expression for individual resistances 
Ri(t) for all moments of time, which reproduces exactly 
the asymptotic behavior 0 for long times and the first 
two terms of ([^ for short times for the total resistance 
R{t). Thus, an approximated expression for Ri{t) can be 
obtained substituting Eq. ([^ into Eq. Assuming 

that 6V <C V, which is the optimal condition for the 
synchronization, one can obtain 


Ri{t) — Roff 


PhSV 

N 


^ Roff 



SV 


+Roff 



SV - 


PiV 

R„//« 


t 


( 12 ) 


The exponential (last) term in the second line of Eq. 0 
decays on a short time scale. Clearly, the ratio of this 
short time scale to the total switching time T, Eq. 0, 
is SVjV < 1. 

Moreover, it is easy to notice that the first two terms 
in the right-hand side of Eq. (12) are dominant at long 
times. These terms are the same for all memristive sys¬ 
tems (they do not depend on the system index i) and 
therefore the memristances of the memristive systems in 
the network are nearly the same. This observation con¬ 
firms our numerical results (see. Fig. [^b)). 

Conclusion.- In conclusion, we have discovered an 
interesting synchronization effect taking place in one¬ 
dimensional memristive networks with elements charac¬ 
terized by a distribution of switching constants. When 
the switching occurs from the high to low resistance state, 
the systems with larger switching constants slow down 
their switching as the voltage falls across these systems 
decrease faster compared to voltages across the systems 
with lower switching constants. As a result, the switch¬ 
ing of all memristive systems occurs coherently with the 
same effective rate regardless the specific switching con¬ 
stants of individual systems. This simple picture explains 
the mechanism of the synchronization effect that is most 
pronounced when the applied voltage slightly exceeds the 
combined threshold voltage of memristive systems. We 
have also demonstrated that the network switching time 
is independent on the number of memristive systems (for 
an appropriately scaled applied voltage) and is defined 
by the harmonic mean of switching constants. 

This work has been supported by the NSF grant No. 
ECCS-1202383. 
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